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From a Monte Carlo study of the ferromagnetic Kondo lattice model for doped manganites, including 
the antiferromagnetic superexchange interaction (Jaf), we found that the ferromagnetic ordering 
was suppressed as Jaf increased. The ferromagnetic transition temperature T c , as obtained from a 
mean field fit to the calculated susceptibilities, was found to decrease monotonically with increasing 
Jaf- Further, the suppression in T c scales with the bandwidth narrowing induced by the antiferro- 
magnetic frustration originating from Jaf- From these results, we propose that the change in the 
superexchange interaction strength between the tig electrons of the Mn ions is one of the mechanisms 
responsible for the suppression in T c observed in manganites of the type (Lao.7- !/ Prj / )Cao.3Mn03. 



PACS numbers: 75.30.-m, 75.30.et, 71.10.-w 
I. INTRODUCTION 

Since the discovery of colossal magnetoresistance 
(CMR) in the doped manganites, Lni— x A x MxLOz 
(Ln=trivalent lanthanide, ^4=divalent alkaline), great ef- 
forts have been made to understand the metal-insulator 
and/or ferro-to-antiferromagnetic transitions driven by 
the substitution of cations, with different sizes, such 
as (La .7_yPrj / )Cao. 3 Mn03.lll _ B For the La . 7 Cao.3Mn0 3 
system, i.e., y — 0, which corresponds to a hole doping 
of x ~ 1/3, a sharp resistivity peak is observed near the 
transition temperature T c « 260 K. While the compound 
is insulating and paramagnetic above T c , it becomes a 
metallic fercomagnet below T c . The double exchange 
(DE) modela has been widely used to describe the na- 
ture of the metallic and the ferromagnetic (FM) states 
in a small range of doping concentrations near x ~ 1/3. 
However, it was claimed thai other effects such as 
Jahn- Teller (JT) distortionjm the orbital degeneracy, 
and the antiferromagnetic (AF) superexchange interac- 
tion (J^)lTE!1 should be included for a more accurate 
description over a wide-range of x values (0 % x ^ 1). 

Recent experimentsETcl on L^/aA^MnOa showed 
that the substitution of trivalent ions such as Y, Pr, and 
Dy at the La 3+ sites having a smaller average radius, 
(Ro), for the lanthanide led to the modification of the 
Mn-O-Mn bond angle (6) and the Mn-0 bond length. 
Those experiments also revealed that these effects caused 
a decrease in T c and a reduction in the magnetization be- 
low T c , as well as an increase in the CMR effects near T c . 
In general, the narrowing of the electronic bandwidth is 



considered to be the origin of the T c suppression effects. 
In this regard, Millis et alB and Rodriguez-Martinez and 
AttfieldQ suggested that the electron-phonon coupling 
arising from the dynamic JT distortion due to the varying 
(Ro) could induce the bandwidth narrowing and suppress 
T c rapidly. However, it was suggested that mechanisms 
other than the elect ron^phjonon coupling could also ex* 
plain the same behavior nu In particular, Radaclli et al.u 
observed that the suppression of T c was more sensitive to 
the bandwidth narrowing due to the change in rather 
than to the electron-phonon coupling. The change in the 
bond angle will not only influence the DE hopping in- 
tegral (t° E ) between e g (Mn)-2p CT (0)-e 5 (Mn) orbitals of 
a FM nature but also affect the superexchange hopping 
integral (t^ F ) through a hybridization of Ihc t2 g (Mxi)- 
2p 7r (0)-t2 g {Mn) orbitals of an AF natures Thus, it is 
expected that both the double exchange and the superex- 
change terms will compete with each other and play a 
significant role when structural distortions occur. There- 
fore, it would be of great interest to study the role of Jaf 
in the doped perovskite Mn-oxides. 

In this paper, we present results of our Monte Carlo 
calculations for the FM Kondo lattice model which was 
modified by introducing the Jaf- From the calculations 
of the spin correlations, we found that the ordered FM 
state changes into an incommensurate (IC) state due 
to the magnetic fluctuations induced by the increase of 
Jaf- It was also shown from the obtained susceptibili- 
ties that the mean-field T c monotonically decreased with 
Jaf- From the density-of-states calculations, we were 
able to determine that the bandwidth narrowing induced 



by the AF frustration originating from Jaf scaled with 
the reduction of T c . This behavior is consistent with 
that seen during recent experimental observations of the 
La 2 /3-2/Prj/Ca 1/3 Mn03 system. 



II. MODEL HAMILTONIAN AND 
CALCULATION METHOD 



effects is observed in Lai-^^l^MnOs with x « 1/3 and 
because dramatic FM-AF and metal-insulator transitions 
occur when La ions are replaced by Pr. With increasing 
Jaf, the system runs into a phase separation (PS) region 
where both AF and IC phases coexist. It should be noted 
that the phase separation induced by Jaf is similar to 
thcjifscontinuity observed in the density as a function of 
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The model Hamiltonian studied in this paper can be 
written as a sum of three terms: 
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where the Si represents the localized spin of the ti g elec- 
trons at the site Ri, the ci a is a destruction operator for 
one species of e g fermions. The ty is the nearest-neighbor 
hopping amplitude for the e g electrons, and Jh > is 
the FM Hund's rule coupling between the e g conduction 
electrons and the ^g localized electrons. The last term 
represents the AF superexchange interaction between the 
t2 g electrons of neighboring Mn sites. In our calcula- 
tion, we treat the localized t2 g spin Si as a classical spin; 
St = 3/2(sin 6i cos fax + sin 9+ sin fay + cos 6iZ). Within 
this approximation, the trace over the e g electrons in 
the partition function can be carried out exactly since 
Eq. (1) is quadratic in the fermionic field. The integra- 
tion over the localized spins {Oi,dirkis performed using 
a standard Metropolis algorithm.EHl3 In our calculation, 
unless stated otherwise, we take t — 1 as the unit of en- 
ergy, Jh — 8t, the temperature T — O.Olt, and L = 6 x 6 
with a periodic boundary condition (PBC) . Because the 
Jaf is estimatedLrEJ to be on the order of Jaf ^ 0.01 
eV while typical values for t and Jh are ~ 0.1 eV and 1 
eV, respectively. The range of values for Jh considered 
in this work looks reasonable in comparison with those 
seen in experiments. 



III. RESULTS 

A. Effects of Jaf in spin-spin correlation 

For a given chemical potential u, = —7.4, the elec- 
tron density (n) (= 1 — a;) was determined as a result 
of simulations. In Fig. 1(a), (n) is plotted as a function 
of Jaf- While (n) stays almost constant for the range 
of Jaf ^ 0.11, an abrupt jump is observed in (n) near 
Jaf ~ 0.115. For small Jaf, the average density is about 
0.65, approximately x ~ 1/3, which is almost indepen- 
dent of Jaf- In this work, we focus on the system with 
x w 1/3 because the FM metallic phase with large CMR 
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FIG. 1. (a) The carrier density vs the antiferromag- 
netic superexchange interaction Jaf on a 6x6 lattice with 
Jh =8t,T- t/100, and (i = -7.4. The notation PS means 
phase separation, (b) Peak intensities S(q) for qfm = (0,0), 
Qic — (cf7T, 7r) (0 < a < 1) and qaf = (tt, n) as a function of 
doping Jaf- 

The spin structure factor 5(q) given by the Fourier 
transform of the spin-spin correlation function, 5(q) = 

Sii eSq ' 1 ('^i ' ^i+i)j was a l so evaluated as a function of 
Jaf and is plotted in Fig. 1(b). The different q values of 
(0, 0), (aw, w) with (0 < a < 1), and (w, w) correspond to 
the FM, IC and AF phases, respectively. The FM fluc- 
tuation is dominant for Jaf < 0.045 and is drastically 
reduced for Jaf ^ 0.05, where the IC component starts 
to develop. Strong IC spin correlations become signifi- 
cantly pronounced with a further increase of Jaf- The 
presence of IC state arising from the frustration of AF in- 
sulating state has been pointed out to play an inmortant 
role in the transport properties in manganites.E3 Also, 
near Jaf ~ 0.115, a small region of Jaf exists where the 
IC and the AF fluctuations are competing. Due to the 
PS effect mentioned above, the two phases should seg- 
regate, leading to the IC and the AF phases. However, 
to understand the fine details of the magnetic proper- 
ties, the one-orbital model is not sufficient due to the 
highly anisotropic nature of the orbital degree of free- 
dom in undoped manganites. Certainly a two orbital 



model is needed to describe the non-trivial interplay be- 
tween the sain and orbital ordering in insulating phase 
of LaMn0 3 fl 
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FIG. 2. Density-of-states JV(w) as a function of e — /i at 
three different temperatures with Jaf = 0.0 at x ~ I /3 with 
a PBC. In the lower panel, the result for a larger system size 
of L = 12 x 12 is shown. 



B. Bandwidth narrowing 

As shown in Fig. 2, we calculated the density-of-states 
N(u>), from the eigenvalues of the Hamiltonian in Eq. (1) 
for many configurations at three different temperatures 
and system sizes with Jaf — 0.0. The broken line at 
u> — fi = indicates the Fermi energy Ept3 It can be 
seen that the bands are split into two major parts with 
the separation between the lower and the upper bands 
corresponding to the energy splitting 2 Jh by the Hund 
coupling. In the low-temperature limit (T = 0.01), the 
spiky shape of each band is due to the finite-size effect, as 
shown by the calculations for two system sizes, L = 6 x 6 
in the upper panel and 12 x 12 in the lower panel. As the 
temperature increases, this shape is smeared out, and at 
the same time the width of each band is reduced due to 
the disorder of the local t^g spins. The system is FM 
metallic since the density-of-states at Ep remains finite. 
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FIG. 3. Density-of-states as a function of e 
Jaf = 0.02, 0.05, 0.08 and 0.11. 
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The N(w) of the lower band is shown for various val- 
ues of Jaf m Fig- 3. These curves are smoothed out 
with a Gaussian broadening of Aw = 0.05. This broad 
spectrum of N(u>) is also observed for < Jaf ^S 0.04, 
which corresponds to the FM region in Fig. 1(b). With 
further increase in Jaf , we find that a significant differ- 
ence exists in the nature of each band between the FM 
and the IC states. For 0.05 < Jaf ^ 0.11, together with 
a strong band narrowing, the wide distribution of eigen- 
values merge into three peaks while the Ep lies at the 
dip of N(u). 

In Fig. 4, the bandwidth determined from the density- 
of-states is plotted as a function of Jaf- We define the 
bandwidth of the lower band as AW = (e max — e m , n ) 
where t max and e m j„ are the maximum and the minimum 
energies of each band in Fig. 2. Surprisingly, the band- 
width has a unique dependence on Jaf in each magnetic 
region, as discussed above. In the metallic FM region, 
the bandwidth is independent of Jaf and is equal to 
about 8. This values is close to the one from the disper- 
sion relation in the 2-dimensional DE model. The most 
dramatic change of W with Jaf is found in the metallic 
IC region. W monotonically decreases by half as Jaf in- 
creases from 0.04 to 0.11. The evolution of the bandwidth 
depicted in Fig. 4 is of great importance since the sup- 
pression of T c in the doped manganites is closely related 
to the bandwidth narrowing. It implies that the AF frus- 
tration arising from Jaf can be an another candidate for 
the band narrowing mechanism in the_doped manganite 
systems besides the JT-polaron effectsO and .the random 
potentials due to the disorder of the cations H 
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FIG. 4. Bandwidth as a function of Jaf for Jh = 8 and 
T = 0.01 with a PBC. 



C. Suppression of ferromagnetic ordering 

The temperature-dependent magnetic susceptibility is 
calculated from X (T) = (S^)/k B TH In Fig. 5, the in- 
verse x(T) is plotted vs temperature for various values of 
Jaf', a Curie- Weiss behavior is seen for x(T) at high tem- 
peratures. The mean-field transition temperature Tm is 
estimated by extrapolating the high-temperature data to 
the low-temperature limit. It is found that Tm decreases 
from positive to negative values, indicating a change of 
the interaction nature from FM to AF, with increasing 
Jaf- The change of Tm with Jaf implies that the mag- 
netic interaction via the DE hopping mechanism is weak- 
ened by increasing the Jaf interaction. 
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FIG. 5. Inverse susceptibility vs temperature for various 
values of Jaf with J H = 8t, T - t/100, and fj, - -7.4. 

Based on the estimated results of Tm in Fig. 5, the 
mean-field temperature vs Jaf phase diagram is shown 
in Fig. 6. Tm in the FM region is determined for two 
boundary conditions; a PBC (circles) and open boundary 



condition (OBC) (squares). The AF region (triangles) 
is determined using only a PBC. The ferromagnetic Tm 
decreases linearly within the error limit, with increasing 
Jaf and reaches zero at Jaf ~ 0.08 above which the AF 
Tm start to develop. The inset shows an experimentally 
proposed phase diagranxl for Lao.7- a Prj,Cao.3Mn03 as a 
function of Pr doping. A decrease in T c when Pr doping 
increases implies a weaker DE hopping overlapping and 
a narrower bandwidth. The two phase diagrams shown 
in Fig. 6 exhibit a similarity between the Pr doping level 
and Jaf- As pointed out earlier, the substitution of Pr 
at La sites causes a smaller (R ) and a smaller bond an- 
gle 0, and the change in modifies both the t® E due to 
the nature of dpa hybridization and the t^ F due to dpir- 
type hybridization. However, while t is more sensitive 
to a change of bond angle 9 due to the nature of dpa 
hybridization, the superexchange term Jaf, which re- 
lies on a <ip7r-type hybridization between t2g(Mn)-p Tr (0)- 
t2 g (Mn), is relatively insensitive to a change in the bond 
angle. Thus, an increase in Jaf, relative to the DE hop- 
ping integral t, is expected when decreases. This be- 
havior is consistent with the experimental evidencefl seen 
in the reduction of the magnetization with (R ). 
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FIG. 6. Phase diagram of the 2D FM Kondo lattice model 
with Jh = St, as obtained by measuring the inverse suscepti- 
bilities at x ~ 1/3. The circles (squares) and triangles repre- 
sent the phase boundaries between the paramagnetic (P) and 
the FM phases with a PBC (OBC) and between the P and 
the AF phases with a PBC respectively. The inset shows the 
results for Lao.7-j/Pr y Cao.3Mn03 from Ref. 1. The solid and 
the dash lines are guides for the eyes. 

Indeed, when the bandwidth in Fig. 4 is com- 
pared with T c in Fig. 6, the suppression in T c for 
Lao.7_j,Prj,Cao.3Mn03 correlates well, at least qualita- 
tively, with the narrowing of the bandwidth. Further, 
the absence of internal spin structure, i.e., IC ordering, 
in Fig. 6 for Jaf ^ 0.08 is due to the mean-field approx- 
imation in determining Tm for X- Thus, the bandwidth 



narrowing in the metallic IC state is thought to be re- 
sponsible for the observed T c reduction. Consequently, 
in addition to the electron-phonon interactions arising 
from Jahn-Tcll distortion, the magnetic frustration aris- 
ing from the enhancement of Jaf should be considered 
as a cause of the bandwidth narrowing leading to the 
suppression in T c . Here, it is noted that the effect of or- 
bital degeneracy is not taken into account in this work. 
Since, however, many optimally doped metallic systems 
near x ss 0.3 are pseudo-cubic in structure, the conduc- 
tion electron hopping can be considered to be isotropic. 
On the other hand, the non-trivial interplay between spin 
and orbital degree of freedom in connection with orbital 
degeneracy has been pointed out to play an important 
role in the spin and orbitakordering in the undoped in- 
sulating phase of LaMn03.Q 



IV. CONCLUSION 

In conclusion, it is suggested that the superexchange 
interaction in the modified FM Kondo lattice model is 
one of the mechanisms responsible for the observed T c 
reduction in Lao.y-yPryCao.aMnOa. With a change in 
Jaf, magnetic frustration is induced and leads to a more 
complicated magnetic phase diagram. The bandwidth 
narrowing occurs for 0.05 < Jaf ^ 0.11 in the metallic 
IC state and is thought to be responsible for the T c re- 
duction in hao^-yPiyCao^MiiOz. With our model, we 
found that bandwidth narrowing, accompanied by T c re- 
duction, can be induced by Jaf without electron-phonon 
coupling. 
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